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This paper proposes a new volume function for calculation 
of the entropy of planar granular assemblies. This function 
is extracted from the antisymmetric part of a new geometric 
tensor and is rigorously additive when summed over grains. 
It leads to the identification of a conveniently small phase 
space. The utility of the volume function is demonstrated 
on several case studies, for which we calculate explicitly the 
mean volume and the volume fluctuations. 
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It has been shown experimentally and computationally 
[1] [2] that jammed granular systems can be described by 
statistical mechanics in appropriate circumstances, vali- 
dating theoretical concepts [3] [4] . The simplest approach 
[3] involves the introduction of compactivity, X — dV/dS 
which plays the role of temperature in thermal systems. 
To quantify X one needs to calculate the entropy S as a 
function of the volume V and therefore the volume as a 
function of the position and coordination of the N grains, 
i.e., a function W to complete the analogy between ther- 
modynamics of equilibrium and these non-equilibrium 
athermal systems E -> V ; H -> W ; S(E, V, N) -> 
S(V,N). This Letter follows a recent analysis [6] of pla- 
nar assemblies in terms of loops and voids. Each grain, 
g can be characterised in terms of the Z g grains which 
it is in contact with, g' , and the position of these con- 
tact points, fggi (see figure 1). For each grain we define 

a centre, f g — (1/Z g ) 5Jg'=i ?gg'i an d vectors fi g that 
connect the contact points. The latter vectors form a 
loop around grain g that is defined to circulate in the 
clockwise direction. Each vector along this loop can be 
uniquely identified in terms of the grain g and a neigh- 
bouring void I. The vectors ri g also form polygons around 
the voids, whose edges circulate in the anticlockwise di- 
rection (see figure 1). To each void we assign a centre 
ft = (V^iIEh^ji where the sum runs over the Z\ 
grains that surround void I. Finally, we define a set of 
vectors Ri g = fi — f g that extend from the centre of grain 
g to the centre of a neighbouring void I. This network is 
self-dual to the f network so that for each R vector there 
is an f vector that intersects it. 

We can now define a one-grain geometric tensor C g 

i 

where the sum runs over all the voids surrounding grain g 
and i, j = x,y index Cartesian components. Each term in 
this expression involves only one self-dual pair of vectors 




FIG. 1. The geometric construction around grain g. The 
vectors r; 9 connect contact points clockwise around each grain 
g and give rise to anticlockwise loops I around each void. The 
vectors Ri g connect from grain centres to loop centres. A 
one-grain geometric tensor is defined as C g = Ri g fi g . 

and has a straightforward geometrical interpretation: Its 
antisymmetric part is exactly A g e, where e = (®i ) and 
Ai g = irig Ri g cos cei g is the volume of the quadrilateral 
formed by n g and Ri g , shown shaded in figure 1. Its 
symmetric part measures the deviation of this quadrilat- 
eral from a perfect rhombus [6]. Note that this holds 
even when the grain has only two contacts, in which case 
the quadrilaterals degenerate into triangles. In an iso- 
static assembly of rough grains (Z g ) = 3, giving rise to 
3JV quadrilaterals altogether. The volume of the entire 
system can then be written as a sum over all grains, 

g 

where / is the identity matrix. It can also be recast 
more conveniently as a sum over all the quadrilaterals 
(henceforth indexed by n) that the f—R pairs make, 

W = -^TnRn cos a n (3) 

n 

where r n — \ fi g \. This volume function (VF) is additive 
as required and its simple form makes it convenient for 
analytical calculations. 

The partition function can be written as [3] 

Z = f e-^efc})!!^ ( 4 ) 
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where {q n } is a set of internal degrees of freedom, 
0({g„}) is a probability density function (PDF) that is 
subject to the contraint that appropriate grains are in 
contact. The coefficient A is the analogue of Boltzmann's 
constant with XX {=1/(3 henceforth) having dimensions 
of volume (L 2 in two dimensions) . The analogue of free 
energy is the effective volume Y = — In Z/(3, the (dimen- 
sionless) volumetric entropy is S = (3 2 dYjd(3 and the 
mean volume is (V") = Y + S/ (3. 

Substituting the VF (3) into eq. (4) we obtain for the 
partition function 

_ /*oo roc plt/l 

Z = T\ l dr n I dRn / da n 

n \_J0 JO J-tt/2 

x e-^ r " R " cosa ^ 2 Q({r n },{R n },{a n }) (5) 

The values of the variables R, r and a are constrainewd 
by a minimal volume that the assembly can attain, v m i n , 
below which further compactification is impossible with- 
out causing grains to overlap, and a maximal volume, 
v m ax, above which dilatation is impossible without the 
loss of mechanical equilibrium and fluidisation of the as- 
sembly. Although these variables are correlated (see be- 
low) it is instructive to first assume independent distri- 
butions. 

A. Correlation-free systems: When the degrees of 
freedom are independent, O can be written as a product 
of probability density functions (PDFs) of the individual 
variables. As a first example, let us consider the follow- 
ing: 



q _ S(Rn - R )S(j n - 70 ) 

n=l 



< r < r n 
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Although simplified, this form yields the essential be- 
haviour observed in experiments [1] [2]. Substituting (6) 
in (5), the one-grain partition function is found to be 



Z l/N = 



-0v n 
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where Z g is the coordination number of the grain and 
the volumes v m i ni v max , and Av are, respectively, 
'mm A i?o7o ' max 

/2, andi?o7o(r max Tmin) /2. The 
volumetric entropy can be directly calculated from this 
expression and the mean volume is 



(v) = Z g 



+ v. n 



1 Av , 
+ - - Ctgh 
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and the volume fluctuations are 



(Sv 2 ) 



1 



AV 



sinh((3AV) 



(8) 
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When [3 — > oo(AA — > 0) we find that (v) — > Z g v m in 
and (5v 2 ) — > 0, while when (3 — > 0(AA — > oo) («) — > 
^ g [(umin + tw)/2 - /?A« 2 /3] and («Ju 2 ) -» Av 2 /3. It 
is interesting to note that expression (8) has the same 
form as the one-dimensional result [3] [5]. The current 
formalism extends that result to two dimensions and 
makes it possible to associate it with a particular dis- 
tribution of degrees of freedom. 

Let us consider a more realistic case 



(Z g )N 
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where <r 7 <C 1, C 7 is the normalisation constant of the 
PDF of 7, and P{r n ) is kept arbitrary for the moment. 
The Gaussian PDF for 7 around 70, whose value is close 
to 1, represents a narrow distribution of the angle a n 
around zero and therefore small devations from a rhom- 
bus [7]. As we argue below, P(R n ) is Gaussian-like and 
is narrower than P(r n ), which justifies its approximation 
as a delta-function. Integration over all 7„ and R n gives 



Z = 



I 



where 

ip{r) = 



drP(r)ip(r)e 
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<j> is the error function and u = (a 7 /jo)(3v. For ucr 7 <C 1 
(high compactivity) 



/ P(r)e 
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where 

x = i- 
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For example, assuming that P(r) is uniform gives that 
the mean volume per grain is 



(V) « z gX 



J max \ ^min 
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and the volume fluctuations are 

(Sv 2 ) = Z, 



1 / XAV V 
{3 2 ^sinh(x/3Ay) J 



(13) 



(14) 



The analysis of these expressions as (3 — * 0, 00 is straight- 
forward and gives the same results as for the previous 
case with Aw replaced by xAw. Similarly, it is straight- 
forwartd to evaluate these quantities for low compactiv- 
ity. 
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The expressions for the volume fluctuations can be now 
related to response functions and diffusion processes in 
granular systems. It is unclear at this stage how large is 
the error introduced by the assumption of independence 
of the variables. One advantage of the VF proposed here 
is that it makes it possible to elucidate this issue. 

B. Correlated systems: To go beyond the assumption 
of independence, the key observation is that correlations 
arise from two sources: (1) relations between the vec- 
tors r and (2) the self-duality of the r- and i?-networks. 
The former can be traced to loops. Whenever M fi g 
vectors close a loop their sum vanishes and one of them 
(2 degrees of freedom) can be expressed in terms of the 
other M — 1 vectors. This introduces two ((-functions 
into the Jacobian, reducing the phase space dimension- 
ality by two. The smallest such loops occur around the 
grains, J2i=i — 0- When these are independent each 
such loop contributed to the Jacobian a term of the form 
Y\i=i P(^i)- The correlations force a modification of this 

term to [Uf^ 1 P(n)]S(f n + E&i" 1 ^)- Note that the 
value of Z g is distributed throughout the system, giving 
rise to a corresponding distribution of such terms in the 
partition function. Voids are also surrounded by loops of 
r vectors (e.g., void / in figure 1), each of which further 
reducing the phase space dimensionality by two. Grain 
and void loops are the only irreducible loops in the system 
(namely, all other loops can be decomposed into these) 
and therefore only they give rise to correlations of the 
first type. This observation has an interesting implica- 
tion: An isostatic systems of TV grains has 3N r-vectors. 
The N granular loops and the N/2 void loops yield 3./V/2 
of these dependent on the rest. Namely, only half the r 
vectors are independent, giving 3N degrees of freedom. 

Turning to the second type of correlations, recall that 
Rig extends from the centroid of grain g, -J- Y^v=i(^g ~ 

l')fi' g , to that of loop I, -^r ^2g'Zi(Zl~ 9')^lg'- Therefore, 
it can be expressed as a linear combination of the vectors 
forming the g and I loops, Ri g = ^felt^ 9 2 a nr n , and is 
uniquely defined in terms of the r-network. On average, 
an R vector depends on (Z g + Zi — 2) = 7 r vectors and, 
since there are two grains to a void, then (R) w V2(r}. 
It follows that P(R) can be safely approximated by a 
Gaussian around this value. The dependence of the R 
variables on the r's further means that the 7 variables 
can also be expressed in terms of these because j n — 



1 — R n ■ r n . In fact, 7„ depends on average on eight r 
variables and therefore P{"fn) can also be approximated 
as a Gaussian. Combining all the above, the VF can be 
written as 



W 



3N/2 

£ ■ 

n,m— 1 



(15) 



The coefficients < a nm < 1 are rational and form a 
sparse matrix of zeroes on the diagonal and, on average, 



seven finite elements in each row. Although appealing, 
this quadratic form is only useful with particular forms 
of P{{r}) and even then it requires knowledge of the 
statistics of a nm . Its main disadvantage is in mixing the 
basic quadrilateral volume units and so losing the clear 
geometrical interpretation of eq. (3). 

Although the correlations reduce the phase space di- 
mensionality to 3N, it is difficult to take them all explic- 
itly into consideration. To make progress, we consider 
all the original 9N variables but include the lowest order 
correlations - those coming from the intragranular loops, 
which are the smallest. The justification for this approx- 
imation is that since the R and 7 variables depend on 
several r degrees of freedom then they are more narrowly 
distributed and so can be considered as background fluc- 
tuations. We believe that this approximation captures 
the correct physics in granular assemblies. 

To illustrate the effects of the intragranular loops, con- 
sider a triply coordinated circular grain of diameter D. 
The three contacts on the circumference form a trian- 
gle of vectors fi g and we wish to know the probabil- 
ity P{r\,r 2 ,r 3 ) that the lengths of its sides lie inside 
(ri,ri +dn), (r 2 , r 2 + dr 2 ), and {r 3 ,r 3 +dr 3 ). The PDF 
of one of the sides falling between r and r + dr is 



P(x = r/D) = 1/ (n^l-x 2 ^ 



(16) 



The first two sides can be chosen independently, but once 
these are in place, the third side is determined by 



= f(x 1 ,x 2 ) = 2 Xl x 2y J '(1 - x\)(l - x\) 

+ {xl + xl)-2x\xl (17) 



and therefore 



P{x u x 2 ,x 3 ) 



S(x 3 - f( Xl ,x 2 )) 

ttV(i -*?)(!-*!) 



(18) 



Although this analysis can be extended to Z g > 3 and 
to non-circular grains, we do not give it here. We next 
demonstrate the effect of intragranular correlations on 
the effective volume. 

Consider a large ensemble of N randomly assembled 
monodisperse circular grains, each contacting exactly 
three neighbours. Ignoring void loop correlations makes 
the grains effectively independent but not the quadri- 
laterls. We take a Gaussian form for P(-f) and practically 
a delta-function for R. The one-grain partition function 
(z = Z X ' N ) is then 



:= / || < "< 2 P(R n )P(-f n )dr n dR n dj n ] 



= [ n [ e-/:r "' t: 

■* n=l 

P(ri)P(r 2 )S(r 3 -f(r 1 ,r 2 )) (19) 

where the lengths r n are measured in units of the grain 
diameter D. Integrating over the R and 7 variables yields 
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Compacitivity X 
FIG. 2. The volume (left axis) and density (right axis) of 
a trivalent granular assembly of monodisperse circular grains. 
The moduli of the intragranular vectors connecting the con- 
tact points are taken to be correlated in the calculation. The 
PDFs of 7 and R are, respectively, Gaussian and a 5-function 
at the periodic lattice value Ro = 2D/\/3. Two density 
plots are presented: one assuming independent r variables 
(squares) and the other taking the intragranular correlations 
into account (+). For X > Xq — 0.265 the correlations com- 
pactify the system while below this value they lead to an op- 
posite effect for the reason discussed in the text. A significant 
feature is that the difference between the two plots is quite 
small. The densities of the two calculations differ only by 
0.009 (« 2%) at AX = 5 and by 0.021 at XX = 0.04 (« 4%). 



dridr 2 P{r 1 )P{r2)ip{r 1 )ip{r 2 )ip{f{ri, r 2 )) x 



exp 



r f-Bl^jrl+rl+f 2 )) _ l3R 0l0 { ri +r 2 +f) 



(20) 



Substituting from (18), we can now compute Y and from 
it the mean volume and its fluctuations. The mean vol- 
ume is plotted in figure 2 as a function of XX. The 
asymptotic value as X — * oo is the average of the close 
and loose random packings, while the value at X — ► 
is the close random packing volume. In terms of grain 
number density, these are 0.473 and 0.537, respectively. 

We can now assess the assumption of independence. 
Using form (18) for P(r) of all three quadrilaterlas, we 
also plot in figure 2 the density of the system with and 
without correlations. It can be observed that the correla- 
tions reduce the density of the system at high compactiv- 
ity, but increase it compared to the correlation-free sys- 
tem at low compactivity. The crossover is at X — 0.265. 
This is expected because the sides cannot all assume si- 
multaneously either r max or r m j n . This constrains the 
high compactivity limit from overexpanding and the low 
compactivity regime from overshrinking. In spite of this 
fundamental difference, it is rather suprising how close 
the densities of the two systems are over the entire range. 
The densities of the correlated and uncorrelated systems 
in the high compactivity limit are, respectively, 0.473 and 
0.464 (1.9% difference), while at the low compactivity 
limit (AX = 0.04) the values are 0.537 and 0.558 (3.8% 
difference). This suggests that simple models that do not 



take the correlations into account may prove sufficiently 
accurate for various purposes. 

To conclude, we have proposed a new VF for planar 
granular assemblies that is exactly additive over grains. 
Nevertheless, we have identified the quadrilaterals, rather 
than the grains, as the basic units of the system because 
it is their volumes that cover the entire volume. A sum 
over grain volumes without involving the quadrilaterals 
has the disadvantage of requiring an assumption on ap- 
proximate grain volume. We point out that the VF pro- 
posed here applies to any assembly in mechanical equi- 
librium, not only to isostatic ones. We have found that 
the relevant phase space is of relatively small dimension- 
ality and identified the sources of geometrical correla- 
tions in the system. We have calculated the effective vol- 
ume and volume fluctuations of several model assemblies 
both with and without correlations. The identification of 
quadrilaterlas as the basic building blocks is relevant to 
modelling. For example, a two-volume model for grains, 
Vi and V2 with occurrence probabilities pi and p 2 , gives 
a mean volume of [5] 



N 



Vl 



— - + Avtanh ( 0Av + ln, — 

2 V V P2 . 



(21) 



Applying the same model to quadrilaterals gives a similar 
expression, but it corresponds to a wider distribution of 
grain volumes both because a grain has Z g quadrilaterals 
and because Z g itself is distributed. Thus, quadrilaterals 
as basic units yield more realistic models without invest- 
ment of extra effort. 
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